NACA TN 2145 


NATIONAL ADVISORY COMMITTEE 
FOR AERONAUTICS 


TECHNICAL NOTE 2145 


LIFT-CANCELLATION TECHNIQUE IN LINEARIZED 
SUPERSONIC -WING THEORY 
By Harold Mirels 

Lewis Flight Propulsion Laboratory 
Cleveland, Ohio 


NACA 


Washington 

August 1950 


lU H l 11 




i 


5 tOEHOO 






TECH LIBRARY KAFB, NM 

111 


IB 


NATIONAL ADVTSORY COMMITTEE FOR AERONAUTICS 


TECHNICAL NOTE 2145 


LIFT- CANCELLATION TECHNIQUE IN LINEARIZED 
SUPERSONIC -WING THEORY 
By Harold Mirels 


SUMMARY 

A lift- cancellation technique is presented for determining 
load distributions on thin wings at supersonic speeds. The tech- 
nique retains certain features of the method recently introduced 
by Theodore R. Goodman, while simplifying and generalizing others. 

A' general 'expression is derived for the load distribution over 
a cancellation wing. This expression permits tbe determination of 
lift distributions on wings that cannot be solved by cancellation 
techniques based on tbe superposition of conical flows. The bound- 
ary conditions for either a subsonic leading edge or a subsonic 
trailing edge can be satisfied. Applications of the expression to 
swept wings having curvilinear plan forms and to wings having 
reentrant side edges are indicated. 


INTRODUCTION 

The method of lift cancellation for obtaining the lift dis- 
tribution on thin wings at supersonic speeds was first suggested 
in reference 1. The lift distribution on a given wing is deter- 
mined by canceling excess lift, through the use of a "cancellation 
wing," on a related plan form having a known loading. This approach 
has been applied by several authors (for example, references 2 to 4). 
The expressions provided in reference 1 are applicable for wings 
that can be generated by the superposition of conical fields. 

A procedure is presented in reference 5 for determining lift 
on a more general class of plan forms than can be handled by coni- 
cal superposition. The method utilizes a surface distribution of 
doublets and an inversion by means of Abel's integral equation and 
is equivalent to a lift cancellation. 
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This report, prepared at the HACA Lewis laboratory, retains 
certain features of reference 5 (that is, the use of a surface dis- 
tribution of doublets and an inversion by means of Abel’s integral 
equation), whereas other features are simplified and generalized. 

The simplification consists in eliminating steps in the procedure 
for obtaining lift distributions. The generalization consists in 
determining a solution that can be made to satisfy the boundary 
conditions for either a subsonic leading edge or a subsonic trailing 
edge (Kutta condition) . The method of reference 5 yields only the 
Kutta solution. The lift-cancellation technique developed herein 
is illustrated by several examples. 

In a concurrent investigation (reference 6), source distribu- 
tions and integral-equation formulations have been applied to obtain 
the loading on a special series of cancellation wings. Reference 7 
employs some of these cancellation wings for the determination of 
lift and moments on swept wings. 


THEORY 

The usual assumptions of an inviscid fluid and small perturba- 
tions are made. The velocity field consists of the free- stream 
velocity U (taken in the positive x-direction) plus the perturba- 
tion velocities u, v, and w. The wing boundary conditions are 
specified in the z = 0 plane. 

The local lift coefficient ACp may be expressed in terms 
of Au. That is. 


AC p = 


PB~Pt ; _ 2Au 

q U U 


( 1 ) 


(All symbols used in this report are defined in appendix A.) Inas- 
much as the local lift coefficient is directly proportional to Au, 
Au will be referred to as "lift" in later developments. 


Lift-Cancellation Method 

The lift distribution on a given wing is to be determined by 
canceling excess lift on a related wing with a known loading. The 
method is illustrated in figure 1. The wing for which the lift 
distribution is desired is shown in figure 1(a) The solution can 
be expressed as the two-dimensional wing (fig. 1(b)) minus a can- 
cellation wing (fig. 1(c)). The loading in region I of the 
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cancellation wing equals the loading in the correspondi n g region of 
the two-dimensional wing and the upwash w in region II of the 
cancellation wing is zero. The loading of the two-dimensional wing 
minus that of the cancellation wing satisfies the boundary conditions 
for the flow about the given wing and is the desired solution. 

The fundamental problem in the lift- cancellation method is then 
to determine the lift in region II of a cancellation wing subject 
to the condition w = 0 in this region and with the assumption of 
a known loading in region I. Solution of this problem is presented 
in the following sections. 


Derivation of Lift-Cancellation Equations 

The lift distribution in region H will be expressed in terms 
of quantities in region I. 

Consider the cancellation wing shown in figure 2. The portion 
of the leading edge to the left of the origin coincides with a Mach 
line. The portion of the leading edge to the right 
(designated r = r^(s)) is shown as a supersonic edge, although no 
restrictions as to a subsonic or supersonic edge are imposed. (A 
plan-form edge is subsonic or supersonic depending on whether the 
component of the free stream normal to the edge is subsonic or 
supersonic.) The line designated r * rgts) separates region I 
and region II and is assumed to be subsonically inclined to the 
free stream at all points. This line corresponds to a plan- form 
edge of the wing for which the lift distribution is desired. 

General solution for load distribution on cancellation wing. - 
The upwash field in the z = 0 plane (due to an arbitrary distribu- 
tion of vorticity Au and Av) may be written, from reference 8, 


9 2 

I 

7i 

[(y-y 0 Mv + (x-x q )au] 

^o d ^o 

2 n 


(x-x 0 ) 2 -P 2 (y-y 0 ) 2 

J3?2 


( 2 ) 


The symbol | designates the finite part of an infinite integral, 
as defined in reference 9. Application of the finite-part concept 
to linearized supersonic-wing theory and the evaluation of the finite 
part of an infinite integral are discussed in references 8 and 10. 

For the present, it will suffice to state the fundamental definition 
of the finite part of an integral with a 3/2-power singularity, 
namely. 


i 



4 


MCA TN 2145 



Ax 


[f(z Q )-f(zj]dz, 




3/2 


2 - -^rk (3) 

(*-a) 1/2 


By a transformation to the Mach coordinates of reference 11, 


x ’ I (s+r) r - ft, ( *- p3r) 

y - 5 (s ' r) 8 " Sg i 

elemental area = ^ dr ds 

M 2 J 


(4) 


equation (2) becomes 


M 

v = — — 
8ir 


p 

r\ 

«--») ^ + <™o) g 

L - Q- 

^o ^o 



T 

[(r-r 0 )(s-8 0 [] 3 / 2 


Upon substitution of the limits of integration, as indicated in 
figure 2, 



r 


r\s 

pf N._ 

< 

ii 

?N 



ds Q 

dr 

^o ° , 


V-> 

(s-s 0 ) 3 / 2 

0 

(r-r c )V2 

r l< 8 o> 


n s 


pr 


ds. 


10 


i^ 7 r2 : 


SA<p 


dr„ 


(r-r ^ 



y 


( 6 ) 
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Integrating by parts, noting that' ACp = 0 at r Q = ^(sq), and 
recalling the definition of the finite part (equation (3)) yield 



( 8 ) 

Similarly, reversing the order of integration (vith appropriate 
changes in limits of integration), integrating by parts, and then 
returning to the original order of integration establish the identity 



( 9 ) 
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The right side of equations (8) and (9) are identical. Equation (6) 
can now he written as 


w = - 


M 

8n 


OS 


rvr 


ds. 


o 


JJ2 


ACp dr 0 


ri(s 0 ) 


(r-rj 


372 


( 10 ) 


For points in region U, w = 0 and equation (10) becomes 


0 = 


os 


ds. 


(s-s ) 3 / 2 

0 ° J r 1 (s Q ) 


ACp dr Q 
(r-r ) 3/2 


(11a) 


or 


p s 

G(r,s 0 ) ds Q 

0= J (a-e 0 )3/2 

0 

where 


(lib) 


0(r,s o ) 



(12) 


Equation (lib) is an integral equation for the unknown function 
G(r,s 0 ). The solution (appendix B) is 


0(r,s o ) = 0 


(13) 


Thus, 


0 = 



r 


n r 2( s o) 

r\ r 


tSP dr 0 


AcPidr 0 



r lK> 


, x( r - r o) 3/2 

r l(s Q ) J 

r 2( S o> 


ACp 


IL 


(r-r Q ) 
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or 


r* 


r 2 (s o> 


TTT 75 


r"- 2 ( s 0 ) 


T l< 0 o> 


4^ dr 0 

77775 


(14) 


The right side of equation (14) will he considered known. Equa- 
tion (14) is then an integral equation for A<Pjj. The solution 
(appendix B) is 



(15) 


Equation (15) indicates that the doublet strength in region II, 
namely AP-q, can be obtained by a line integration along s Q = s 

in region I. The geometric interpretation of the various terms in 
equation (15 ) is shown in figure 3. 


It can be shown, by expanding rfPj in a Taylor’s series about 
r = r 2 (s), that equation (15) yields a continuous solution 
(A«Pn = APj) at r = r 2 (s). (A discontinuity in AP implies a 
lifting line (reference 12) and is unrealistic . ) 

The lift distribution in region II can be expressed as 


Au 


II " 


d^II 

5x” 


d Atp H Sr ^ ACp II Ss 

Sr Sx Ss Sx 



or, from equation (15), 


¥ au ii 




pr 2 (s) 

aJ r-r 2 (s) 

Aqpj dr Q 

Jt 

r x (B) (r - r °>4 r 2< s >- r o 


(16) 
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Differentiation yields (See appendix C.) 


Au II - 


aJ r-r 2 (s) 


p r 2 ( s ) 


Au_ dr 
I 2_ 


(s) 


(r-r 0 ) r 2 C s ) " r o 


dr 2 (s) 


ds 


28rt/^r-r 2 (s) 


n r 2 ( 8 ) 

(8Au I -AT I )dr c 

, \ f\j r 2^ s ^~ r o 

T-. (S) * 


(17a) 


Equation (17a) is *the desired expression for the lift distribution 
in region II in terms of quantities in region I. 

Consider Aujj to consist of two components, Aujj 1 and 
AUjj", where Au-^' and Atn." are the first and second terms 

on the right side of equation (17a), respectively. Investigation 
of the integrals indicates that at r = r 2 (s), Aujj* = Au-|-j 

whereas Aujj", in general, has a half-order singularity. 


When region II is to the right of region I (fig. 3(h)), the 
integration for Aujj is conducted along the line r Q = r and 

may he written as 


Aun 


A/ S “ s 2 (r) 

jt 


ps 2 (r) 

Au t ds 

L- o 

(s-s 0 ) /\j B 2 (r)-3 ■ 

s x (r) 


1 - 


ds 2 (r) 


ps 2 (r) 


dr 


2p* Ajs-s 2 (r) 


St (r) 


(pAu I +Av I )ds (; 
/\J s 2 (r)-s 0 


(l7h) 
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Discussion of equations (17a) and (17b). - In the paragraph 
preceding equation (2), the line r = r 2 (s) was described, as sub- 

sonically inclined at all points to the free stream. This condition 
is necessary so that the inner integral in equation (lla) (that is, 
G(r,s Q )) can be equated to zero for all points in region II. If 

this restriction on r = r 2 (s) is not satisfied, the development 
beyond equation (lla) becomes invalid. The derivation of cancella- 
tion equations when r = r 2 (s) is supersonic was not undertaken 

because such problems can be solved more simply by other methods. 

In regard to the boundary conditions, it has been assumed that 
Au-j- is specified. Equations (17a) and (17b), however, indicate 

that a knowledge of Av-j- is also required in order to obtain a 
solution for Aujj. Two possibilities exist, as illustrated in 

figure 4. In the first case (fig. 4(a)), region I is upstream of 
region H (along the line r = r 2 (s)) and Avj is uniquely defined 

by the specified Auj according to the relation 


AVn 


I»r 

“■§7 ' 


J ^(y) 


Aui dx Q 


(18a) 


The integration is conducted along lines of constant y. In the 
second case (fig. 4(b)), region II Is upstream of region I (along 
the line r = r 2 (s)) and the expression, for Avj in region I 

.(for y<0) is 




px 2 (y) • px 

4T I-^ 


Aun dx Q + Au x dx Q 


-Py J x 2(y) 


Equation (18b) indicates that a knowledge of Aujj is required in 
order to find Avj. But Avj must be known (equation (17a)) before 
Aujj can be found. Thus, the solution for Aujj from a specified 
Auj is not unique for the configuration of figure 4(b) and an 
additional boundary condition must be imposed. The line r = r 2 (s), 
however, corresponds to a plan-form edge of the airfoil whose load 
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distribution is desired. The situation indicated in figure 4(h) 
occurs when r = r 2 (s) corresponds to a subsonic trailing edge. 

The additional condition to be imposed is therefore the Kutta con- 
dition. In terms of the cancellation wing, this condition requires 
that the perturbation velocities be continuous in crossing r = r 2 (s). 

Solution for Aujj satisfying Eutta condition at r =r 2 (s) . - 
It will now be shown that when the Kutta condition is imposed at 
r = r 2 (s), the appropriate Avj distribution is such as to make 

the second integral in equation (17a) identically zero; that is. 


O r 2 ( s ) 


r-i(s) 


pAu I -Av I 

h<=k] 1/2 


dr Q = ° 


^ACpp 

y 

r o 


or, inasmuch as PAuj-Avj = M — — , 

pr 2 (s) 




"l(s) 




= o 


This concept and its proof follow from a suggestion of I. S. Kibner 
of the HA.CA Lewis laboratory. 

Thuf, from equation (7), (12), and (13), 

or 


5^° 


r^(s) ( r ~ r o) 


1/2 


= 0 


(19) 


for all points (r,s) in region II. Therefore, 


pr 2 (s) 


"l(s) 


i?: ° 


r 


(r-r 0 ) 


1/2 


Sa«P- 


n 


r 2(s) 


Sr o 

(r-r Q ) 


dr. 


T/2 


( 20 ) 
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Upon taking the limit as r approaches r 2 (s), equation (20) 
becomes 


pr 2 (s) 

r 

r\T 

Sa<P t 


^11 

ST _ un 


-ST 

, . lf2( s )- r o] 1/2 r-^r 2 (s) 

r x (s) *- 


, * ^ 0 ) 

r 2( s ) J 


( 21 ) 


Hovever, 3tfPj/3r 0 must be continuous 'in the vicinity of r 2 (s). 
(The perturbation velocities on the basic wing can be discontinuous 
only along Mach lines op along plan-form edges. Inasmuch as 
r = r 2 (s) is neither of these cases, all derivatives of ACp-j- must 
be continuous in the vicinity of r = r 2 (s).) When the Kutta con- 
dition is imposed, dACpjj/Sr Q is therefore also continuous (and 
bounded) in the neighborhood of r = r 2 (s). Then, using a mean 
value for ^ifpjj/br Q} 


lim 

r-^r 2 (s) 


r r 


^(a) 


~ST ° 

; (r-r 0 )V 2 


lim 

r->r 2 (s) 


— 

r r 

/asPnA 

dr D - 

\ ^ r o/ (r 2 (s)<r 0 <r) 

___ 

(r-ro) 1 / 2 

r 2 (s) 


= 0 


( 22 ) 


Therefore 


pr 2 (s) 


-r — dTr. 

dr„ 0 


r l( .) [ r 2 (6) - r °] V2 


0 


(23) 


which was to be proved. 
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The solution for Au-j-j that satisfies the Kutta co ndi tion at 
r = r 2 (s) is then, from equations (17a) and (23), 


pr 2 (s) 


Aun = 


^r-r 2 (s) 


Auj dr Q 


, , ( r - r o)A/ r 2^)- r c 

for the wing of figure 3(a). Similarly, 

ps 2 (r) 


Au n = 


/^s-Sg(r) 


it 


Au t ds 
I o 


Sn (r) 


(s-s 0 )/\/s 2 (r)-s c 


(24a) 


(24b) 


for the wing of figure 3(b). 

An alternate derivation of equations (24a) and (24b) (appen- 
dix D) indicates that only solutions satisfying the Kutta condition 
will result from the integral equation formulations of reference 5. 

Sidewash in region II. - An expression for Avjj can be obtained 
by differentiating equation (15) with respect to y. The result is 


- V r r 2 


AV H = x 


vJ 


1 + 


dr 2 (s) • 


ds 


2«/^r-r 2 (s) 


pr 2 (s) 

r x (s) 

nr 2 (e) 


AY I 

(r-r D ) Ajr 2 (s)-r 0 


PAUj-Avj 

, l( s) 


dr. 


Si m ilarly, for region II to the right of region I, 
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Av II = 


a/ s-s 2 ( r ) 


ps 2 (r) 


s x (r) 


Av I ds Q 


(s-s 0 ) Ajs 2 (r)-s c 


1 + 


ds 2 (r) 


P s 2^ r ) . 


ds 


2ir^s-s 2 (r) 


s l( r ) 


■ PAU I +AY I ^ 

A/s 2 (r)-s 0 0 


When the Kutta condition applies, these equations become, 
respectively. 


Av II = 


/^r-r 2 (s) 


pr 2 ( s ) 


^(s) 


Av T dr 
_I Q_ 


(r-r D ) /\/r 2 (s)-r c 


and 


^s 2 (r) 

Av I ds o 

Sx(r) (s-s 0 )^s 2 (r)-s 0 

It should be noted that vhen r = rg(s ) corresponds to a subsonic 
trailing edge, Avj, as well as Avjj, is not generally known. 

The preceding expressions are therefore primarily useful for those 
problems where r = r 2 (s) corresponds to a subsonic leading edge. 


Av- 


II 


. /^s-s 2 (r) 


APPLICATIONS 

The loading in region II of a cancellation wing is given by 
the line integrals of equation (17a) or (17b). When the Kutta 
condition is imposed at a subsonic trailing edge, the expressions 
reduce to equations (24a) and (24b). These equations can be used 
to find the load distribution on a large variety of wings that 
cannot be solved by cancellation techniques based on conical 
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superposition. Wings witn curvilinear plan forms or arbitrary 
camber are examples. In each case, however, the solution for the 
related wing must be known. 

The equations are applied in several illustrative examples. 
Only the solution associated with the cancellation ving is con- 
sidered. The complete solution consists in the loading of the 
related wing minus the loading of the cancellation wing. 


Leading-Edge and Side-Edge Cancellations 

In these cases, the lift to be canceled is upstream or to the 
side of the plan form for which the loading is desired (figs, 5 
and 6) . 

Tip region of swept wing. - The loading in the region influ- 
enced by the side edge (II a and 11^ of fig. 5) of a swept wing 
having a subsonic leading and a supersonic trailing edge can be 
obtained by canceling excess lift on a trian g ular wing. The Kutta 
condition is applied across the portion of r = r 2 (s) influencing 

region H^. The lift to be canceled in region I 1 b (reference 13, 
equation (23)) 


, Jgk _ - = » i W> (25) 

Aje^-p 2 ^ ^je 2 ( s+r) 2 -(B-r) 2 

where H and 6 are constants defined in appendix A. The doublet 
distribution in region I a , again from reference 13, is 

A q> = = M/\Je 2 ( s+r) 2 -(s-r) 2 

I a M 

from which 

4V. = SA ^ a = jgk - = - (26) 

a 37 A/aV-p 2 / A/a 2 (^r) 2 -(s-r ) 2 

The sidewash distribution in region 1^ (that is, AVj^) could be 

found by an integration of the type indicated in equation (18b). A 
knowledge of AVj^, however, is unnecessary in the present problem 

because the Kutta condition is applied for region H^. 


o 
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The loading in region H a is obtained by substituting equa- 
tions (25) and (26), vith r replaced by r Q , into equation (17a), 
which yields 


pr 2 (s) 


Au 


n. 


HS^A^r-r^Ts) 


(s+r ) dr 
v o' o 


H 


b - aS 

x ^r-rg(s) 


r l(s) 

pr 2 (s) 

r l(s) 


( r “ r 0 ) A|r 2 (s)-r o ^0 2 (s + r ^)2-(s-r o ) 2 


[j 2 (sfr 0 )+( s -r 0 )]dr 0 

/J r 2 (B)-r 0 /je 2 (s+ro) 2 -(s-r 0 ) 2 


(27) 


For region H^, the Kutta condition applies and 

pr 2 (s) 


Au II t " 


_ KO 2 f^r-r 2 (e) 


n 


(s+r Q ) dr Q 


r l(s) 


(r-r Q ) /\|r 2 (s)-r 0 A/e 2 (s+r 0 ) 2 -(s-r 0 ) 2 

(28) 


Equations (27) and (28) reduce to elliptic integrals of the first, 
second, and third kind upon transforming the variable of integra- 
tion from r Q to ( 0 Q according to the relation 


= TTe [d-oJs + W] 


( 29 ) 


where 


a p = (l+S)r„(s)-(l-e)s 


Equations (27) and (28) may then be written 
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n = — ~ a g = (1+0) r 2 (s)-(l-0)s 

Beentrant aide edge. - A plan form has a reentrant edge if a 
line of constant y intersects the plan form at more than two 
points . 

The load distribution in the region influenced "by. the reentrant 
side edge is to "be determined for the wing of figure 6. The side 
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edge is first, for simplicity, the straight line r = K 2 s, which 
is a subsonic trailing edge across which the Kutta condition is 
applied. The side edge then alternately becomes a subsonic leading 
and a subsonic trailing edge. The load distribution in region I is 

simply the Ackeret value Au-j- = 2s2, and Avj = 0. Regions II a , 

lib, and- He are considered separately. 

Region II a . 

From equation (24a) with Auj = pp 

EgS 

2aU dr Q 

s P(r-r 0 )^K 2 s-r 0 


AU II. 


A/r-K 2 


4ccXJ 

P* 


tan' 




(Kg+I)s 


■- K 2 b 


(32a) 


or, in x,y- coordinates 


Region II t . 



4<xU 

P* 


tan - - 1 - 


x+Py 

p(mgx-y) 


(32b) 


A knowledge of 



is required. From equation (18b), 


= ^ 



^y/mg 


r x 

4dU 

P* 

tan -3 - A 
-Py 

x ° +f V d^ i 2aU 

^(mgXo-y) P 

J 

dxQ 

y/®2 


= -*2L . 2an aIeEI 

^ pmg] V e 2 _1 


( 33 ) 




Region II C . 

Because the Kutta condition is applied, 
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Trailing-Edge Cancellation 

The calculation of lift distributions on swept wings haring subsonic trail ing edges 
requires cancellation wings of the type shown in figure 7. These wings cancel that part 
of the lift of the basic triangular wing that is downstream of the trailing edge of the 
swept wing (references 3, 4, 6, and 7). The lift is specified in region I. The lift in 
regions II and III is to be determined subject to the conditions that w = 0 and that the 
Kutta condition applies at r <= r-^s) and r = r g (s). 


The wing of figure 7(a) differs from the previously discussed cases in that two unknown 
regions (II and III) are continuously interacting. A special treatment is required in order 
to obtain the loading in regions II and III. (See, for example, reference 6.) Approximate 
solutions can be obtained, however, using equations (24a) and (24b). For example, if the 
load at (r,s) in region II is desired, first assume that Aujjj is known. Then 



“ 

nr X ( s) 

pr 2 (s) 


/\|r-r 2 (s) 


! Au t __ dr 

.III o + 

AU] 

. dr 

o 



Mc r (*-r 0 )A/ r 2( s )- r o ^ 

^(s) ( r_r o) A 

r 2 (s)-r 0 


vj 

26 


— 


(36) 


An expression for Aujjj is, by integrating along lines of constant r Q (fig. 7(a)), 


Au 


III 




P s 2< r o) 

P 8 l( r 0 ) 

Ajs-si(r 0 ) 



Auj ds D 

jt 


Mc r ( s - 0 o)^ B l(ro)- B o ^ 

S2 (r 0 ) (S ' So> Ai s l< r o)- s o 


\J 

26 



H 

to 
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An approximate expression for Au^j-j- is then 


An 


III 


- V s " B i tr„) 

M 1 ‘ ** * — 


psi(r 0 ) 


Au T da 
J, S_ 


s 2^ r o) 


(s-s 0 ) A/s 1 (r 0 )-s c 


( 37 ) 


Equation (36 ), which may now be written in terms of Auj by substituting equation (37) 
for AUjjj, becomes 


Au II « 


A/r-r 2 (s) 


r 

n r x (s) 


n B i( r 0 ) 

J 

^0 

/Js- s i(r 0 )' 

AU X ds Q 


Mc r ( r - r o)^ r 2( s )- r o. 

it 

s p (rj ( 8 - s o)^ s l( r o)- 8 o 

V 

"2P 




P r 2(®) 




A u I 




rjU) ' r - r o)4 r 2 (')- : 




(38) 


The first term on the right side of equation (38) approximates the contribution of 
region III to the loading in region II. This term, as indicated in reference 6, is 
negligible for the commonly encountered Au-r distributions (corresponding to steady lift, 

(3 dy-j_ (x) 

roll, or pitch) and > 0.5. For those cases, equation (38) simplifies to 


A63T 


NAjCA TN 2145 



1297 


pr 2 (s) 


Au II ~ 


V r-r 2 (s) 
n 




r ( 8) (r-r 0 )^r 2 ( s )-r c 


' (39) 


The Kutta condition at r = r 2 (s) is satisfied by both equations (38) and (39). Equa- 
tion (39) can be reduced to canonical form by the substitution indicated in equation (29). 


When region I has a partly supersonic leading edge (fig. 7(b)), it is possible to 
■write exact expressions for the linearized load distribution in regions H and HI. For 
example, the load at point (r,s) of figure 7(b) is 


* 


Auji 



r 

pr 1 (s) 


p 8 l( r o> 

Aj r - r 2 ( 8) J 


dr 

o 

A/s-Bi(r o ) 

AUj dB o 

it 

V. 

c (r-r 0 )^r 2 (s)-r 0 

it 

is 3 (r 0 ) 


\r ? (s) 


rx(s) 


Auj dr 0 


(r-r 0 )/\|r 2 (s)-r c 


(40) 


Successive Cancellations 

A cancellation wing may induce lift that itself must be canceled in order to satisfy 
boundary copditionB completely. Thus, in figure 8, the cancellation of lift in region I 
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induces lift in region I*. The cancellation of lift in region I* 
induces lift in region I", and. so forth. Each of these cancella- 
tions is handled, as previously described. These , computations are 
very tedious when lift is induced upstream of a subsonic leading 
edge (for example, region I* of fig. 8), inasmuch as a knowledge 
of the sidewash (Av^,), as well as of the lift distribution 

(Auji), is needed in order to continue the cancellation process. 
Numerical methods are generally required. 

Successive cancellations axe discussed more extensively in 
references 3 and 4. 


SUMMARY OF ANALYSIS AND APPLICATIONS 

A general expression was determined for the lift distribution 
over a cancellation wing. The expression is valid when the plan- 
form boundary (on cancellation wing) separating the region of zero 
upwash from the region for which the lift is specified is everywhere 
subsonically inclined to the free stream. This expression permits 
the determination of lift distributions on wings that cannot be 
solved by cancellation techniques based on conical superposition. 

The boundary conditions for either the flow about a subsonic leading 
edge or a subsonic trailing edge can be satisfied. 

The lift cancellation technique was illustrated for swept wings 
having curvilinear plan forms. Leading-edge and trailing-edge 
cancellations were considered. In addition, the loading in a region 
influenced by a reentrant side edge was found. 


Lewis Flight Propulsion Laboratory, 

National Advisory Committee for Aeronautics, 
Cleveland, Ohio, January 16, 1950. 
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APPEflTDIX A 
SYMBOLS 

0?he following symbols are used in this report: 
(l+0)r-(l-0)s 
(l+0)r 1 (s)-(l-0)s 
(l+0)r 2 (s)-(l-0)s 

PB-Pt 


p 
c 

c r 

E(0,k) 


F^k) 


local lift coefficient, 

H 

constant 

root chord of swept wing 

elliptic integral of second kind, 
pain 0 

a/ 1-k 2 &)_ 2 


E(^,k) = 




elliptic integral of first kind, 
psin <f> 


F(0,k) = 


dd) 


0 l-0) o 2 )(l-t 2 ffl o 2 ) 


G(r,s 0 ) function, of r and s 0 defined hy equation (12) 


H 

K 

k 


2oU 




PE 


slope of plan-form edge in r,s-coordinates, dr/ds 
modulus of elliptic integrals 
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M 


Mach number 

m 


slope of plan- form edge in x,y- coordinates, dy/dx 

n 


parameter of elliptic integral of third kind 

P 


local static pressure 

q 

s= 

IpU 2 

r ' r ol 

s > 8 oJ 

[ 

Mach coordinate system (equation (4)) 

u 


free- stream velocity 

U,T,W 

perturbation velocities in x-, y-, and z-directions 
respectively 

Au 

= 

Urp-Ug (proportional to local lift) 

Av 

= 

V t b 




y,y G 

► 

Cartesian coordinate system 

z > z oJ 



a 


angle of attack 

P 

= 

aJm 2 -i 

8 


semivertex angle of triangular wing 

9 

= 

P tan 8 

n(0,: 

n,k) 

elliptic integral of third kind. 


psin (f> 


au> 


(l+n(0 o 2 ) ^(l-k 2 0) o 2 )(l-0) o 2 ) 
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P 

T 

<P 

ACp 

«o 


density 

area of integration 
amplitude of elliptic integrals 
perturbation velocity potential 
doublet strength, cpij-Cpg 

integration variable 


Regions : 

X 

I a^ I b^ * * * 
II 


region on cancellation wing for which loading is 
specified 

subdivisions of region I 

region on cancellation wing for which w = 0 


II a ,II^, . . . subdivisions of region II 


IH 


additional region on cancellation wing for which 
w = 0 


Special designations: 

r = r^(s) r as function of s along plan- form boundary 1 
s = s-|_(r) s as function of r along plan-form boundary 1 
y = y-j_(z) y as function of x along plan-form boundary 1 
x = x x (y) x as function of y along plan- form boundary 1 
r = r g( s ) r as faction of s along plan-form boundary 2 
s = Sg(r) s as function of r along plan-form boundary 2 


and so forth 
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Subscripts: 

1,2,3 

i ,n 

B 

T 


refers to plan- form boundaries 1, 2, and 3, respectively 
refers to regions I and II, respectively 
bottom surface of z = 0 plane 
top surface of z = 0 plane 


o 


variable of integration 
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APPENDIX B 


SOLUTION OF INTEGRAL EQUATIONS 

Consider the following integral equation (in the notation of 
the appendix in reference 5), where the function f(x) is assumed 
known and the function u( |) is to he determined: 


f(x) 


ox 


u(|) 


ox 


(-5) 3/2 


2u(x) 

( x - a )^/ 2 


(Bl) 


After an Integration hy parts, equation (Bl) may he written 


f(x) 

2 


- jM„ 

(x-a) 1 / 2 


ox 




(B2) 


Equation (B2) is now an integral equation of the Ahel type. The 
continuous solution for u(| ) is (reference 14) 


pz 

u(z) = - % (B3) 

2 ” (z-x ) 1 ' 2 

<J a 

evaluated at z - This result is presented in reference 5. 

Equation (lib) corresponds to equation (Bl) with u(£) =. G(r,s 0 ) 
and f(x) = 0. The solution for G(r,s 0 ), according to equa- 
tion (B3), is then 


G(r,s 0 ) = 0 


(13) 
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Equation (14) corresponds to equation (Bl) with. 
t = r Q a = r2(s c ) 

P r 2( s o) 

x = r f (x) = - 

U 


ri(s c ) 


dr 


u(|) = ACp. 


II 


The solution for A ( pjj according to equation (B3) is then 




ACp ii - s 


P r 2^ s o) 


dr 


t \ 4 ^ 

r 2 (s 0 ) u 


A< ?I ^Q. 

(r-r Q ) 3 / 2 


^i(8 0 ) 

After reversing the order of integration and integrating. 


ACPj-j = 


n T 2(®o) 


A/z-r 2 (s 0 ) 


«Pl ^o 


r l( s o) 


(z-r 0 ) tjr 2 (s 0 )-r c 


Equation (B5), evaluated at z = r and s Q = s, yields 


A<Pii = 


4 r-r 2 (s) 


pr 2 (s) 


APPi dr Q 


*l(s) 


(r-r Q ) A/r 2 (s)-r c 


(B4) 


(B5) 


( 15 ) 


The derivation of equation (15) is similar to that for equa- 
tion (16) of reference 5. 
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APPENDIX C 

DIFFERENTIATION TO OBTAIN Aujj 
T he differentiation indicated in equation (16) 


26 


is to be conducted. 
First 




pr 2 (s) 

(a + 

Ai r- p 2 (s) 

A<p r 

\5r 35/ 

L ” J 

r x (s) ( r " r o)^ r 2( s )- r o 


2p A 
M AU H 


ds 


p r 2 (s) 


2it/Jr-r 2 (s) 


AC Pr ^o- 


n (s) (r - 1 '°W r 2 (S, ' r = 


(16) 


^|r-r 2 (s) / j_ a \ 

pr 2 (,s; 

ACpj- 

it \§r ds/ 

ri(B , (-r 0 )^ 


Inasmuch as A^P-j. is a function of r Q and s. 


"5? 


pr 2 (s) pr 2 (s) 

r l( s) ( ^„)^(.)-r 0 J n(>) 


ACfX n 


(ci) 


(r-r Q ) 2 A/r 2 (s)-r c 


(C2a) 
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and 


f^ 2 ( s ) 


5 

Si 




*»I 


(^K 


ll(s) 


(r-r 0 )^r 2 (s)-r c 


r^e) 


( r-r o) t\J T 2( B )~ T c 


Urn 

T o~> r 2< s) 


1 1 ^ 2 ( 6 ) 

o 

4 

< 

o 

5 L 

l 2 J 

ri(s) (r - r o)[ r 2 ( s )- r o ] 3/2 


( A<P I^r 0 = r 2 (s) dr 2 (s) ( AC Pl)r 0 = r x (s) cLri(s) 

(r- r o )^r 2 (s)-r o 63 | [r-r^ 8 )] ^r g (s )-r x ( s ) 48 


(C2b) 


However, ( AC Pi)r 0 = ^(s) = 0 integration by parts. 


pr 2 (s) 

A<P t dr 

-I— — Q _2 

pr 2 (s) 

1 1 

*~b 

_h ^ 
9- ■ 

<i 

l>° 

i^.L i 

, y ( r - r of r 2( S )- r 3 3/2 

vJr-, (s) L ^ vJ 

*l(s) 

(r-r 0 )^r 2 (s)-r 0 


dr„ 


lim. 


2A<p, 


r o~ )r 2^ s ^ |^(r-r 0 )/^r 2 (s)-r c 
so that equation (C2b) can be written as 


J r l(s) 
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pr 2 (s) 


pr 2 (s) 


5 

Si 




:i(S )( r - ro )V» 


SAqpj- 


sr ^ 


r^.) (r-r 0 )^T 2 (s)-r c 


dr 2 (s) 

pr 2 (s) 

dA^ _ A<Pj- 
_dr Q + (r-r Q )_ 

^o 

ds 

ri(s) 

(r-r 0 )^r 2 (s)-r 0 


( C3 ) 


Upon substituting equations (C2a) and (C3) into equation (Cl) and 
integrating by parts those integrals containing equation (Cl) 

finally reduces to 


t\J r-r 2 (s) 


AU II = 


1 - 


dr 2 (s) 


r 


^ 2 (s) 

ri(s) 
r 2 (s) 


Au t dr 
I o 


( r - r o)^/r 2 (s)-r c 


ds 


2 r_r 2^ s ) 


(PAu t -Av t ) dr 
I I' o 


•l(s) 




(17a) 
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APPENDIX D 


ALTERNATE DERIVATION OF SOLUTION SATISFYING 


KUTTA CONDITION AT r = r 2 (s) 


The integral equation formulation in terms of A<p (equa- 
tion (14)) resulted in a solution that was continuous in Acp 
(equation (15)) but discontinuous, in general, in the derivative 


dAqp 


= Au (equation (17a)) at r = r 2 (s). 


In order to obtain a 


solution continuous in Au, an integral equation may be formulated 
that is similar to equation (14), but in terms of Au rather than 
A<P. The inversion shown in appendix B should result in a solution 
that is continuous in Au but discontinuous in the derivatives of 
Au at r = r 2 (s). 


Consider equation (10) for the w distribution in the 
z = 0 plane. This equation will be differentiated with respect to 
x using a technique introduced in reference 15 (equations (l) to 
(3)). The expression for w at any point (r,s) Is, from 
equation (10), 



where T is the area abc in figure 9(a). The wing is moved 
upstream a distance dx (fig. 9(b)), keeping the coordinate system 
fixed in space. The expression for the upwash at (r,s) now 
becomes 
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The second term on the right side of equation (D2) is zero because 
ACp = 0 along the leading edge. Subtraction of equation (Dl) from 
equation (D2) then yields 


5v JM_ 

^x 8 k . 


Au dr Q ds Q 
(«-•„ ) 3/2 (r-r H 2 


For points in region II of a cancellation wing, 
for the wing of figure 3(a), 


dw 

Si 


(D3) 

0. Thus, 



r % 

nx 

Au dr 

o 

J 


■ f 

r l( s o> 


(D4) 


This equation is the same as equation (11a) except that Au replaces 
A<P. The inversion by Abel's integral equation, for Au-j-j in terms 

of Auj then gives (from equation (15)) 


pr 2 (s) 


AUn = 


hjr-rAs) 




r l(B) ( r ~r 0 ) ffeisT- 


(24a) 


Inasmuch as the integral equations of reference 5 are formulated 
in terms of Au and are inverted by means of Abel's integral equa- 
tion, only solutions satisfying the Kutta condition will result 
therein. 


references- 

1. Lagerstrom, Paco A.: Linearized Supersonic Theory of Conical 

Wings. NACA TN 1685, 1948. 

2. Coleman, T. F.: Supersonic Lift Solutions Obtained by Extending 

the Simple Linearized Conical Flow Theory. Eep. No. CM-440, 
North American Aviation, Inc., Feb. 20, 1948. (Bumblebee Proj. 
Contract NOrd 9784.) 

3. Cohen, Doris: The Theoretical Lift of Flat Swept-Back Wings at 

Supersonic Speeds. NACA TN 1555, 1948. 



34 


NACA nr 2145 


4. Cohen, Doris: Theoretical Loading at Supersonic Speeds of Flat 

Swept -Bach Wings with Interacting Trailing and Leading Edges. 

NACA TN 1991, 1949. 

5. Goodman, Theodore B.: The Lift Distribution on Conical and 

Nonconi cal Flow Begions of Thin Finite Wings in a Supersonic 
Stream. Jour. Aero. Sci., vol. 16, no. 6, June 1949, 
pp. 365-374. 

6. Bibner, Herbert S.: -Some Conical and Quasi-Conical Flows in 

Linearized Supersonic Wing Theory. NACA TN 2147, 1950. 

7. Bibner, Herbert S.: On the Effect of Subsonic Trailing Edges 

on Damping in Boll and Pitch of Thin Sweptback Wings in a 
Supersonic Stream. NACA TN 2146, 1950. 

8. Bobinson, A.: On Source and Vortex Distributions in the Linearized 

Theory of Steady Supersonic Flow. Bep. No. 9, College Aero. 
(Cranfield), Oct. 1947. 

9. Hadamard, Jacques: Lectures on Cauchy's Problem in Linear Partial 

Differential Equations. Oxford Uhiv. Press (London), 1923, 
pp. 133-135. 

10. Heaslet, Max. A., and Lomax, Harvard: The Use of Source-Sink 

and Doublet Distributions Extended to the Solution of Boundary 
Value Problems in Supersonic Flow. NACA Bep. 900, 1948. 
(Formerly NACA TN 1515.) 

11 . Eward, John C . : Use of Source Distributions for Evaluating 

Theoretical Aerodynamics of Thin Finite Wings at Supersonic 
Speeds. NACA Bep. 951, 195° 

12. Mirels, Harold, and Haefeli, Budolph C.: Line-Vortex Theory for 

Calculation of Supersonic Downwash. NACA TN 1925, 1949. 

13. Heaslet, Max A., and Lomax, Harvard: The Calculation of Down- 

wash behind Supersonic Wings with an Application to Triangular 
Plan Forms. NACA TN 1620, 1948. 

14. Whittaker, E. T., and Watson, G. N.: Modern Analysis. The 

MacMillan Co. (New York), 1943, p. 229. 

. Eward, ‘John C.: Theoretical Distribution of Lift on Thin Wings 

at Supersonic Speeds (An Extension) . NACA TN 1585 , 1948 . 


15 



LU 


w=-au 

u=Auj 


/ w=-*aU 
Au=Aut 


w=-aU 

Au=Aut 


(b) Two-dimensional wing 


Au=Au 


(c) Cancellation wing. 


Figure 1. - Superposition to obtain lift on given wing by 
canceling lift on two-dimensional wing. (Given wing equals 
two-dimensional wing minus cancellation wing. ) 
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Figure 8. - Successive cancellations 
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